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Abstract 

Recent data on diffractive photoproduction of dijets are analyzed within the 
framework of fracture functions and paying special attention to the consequences 
of the use of different rapidity gap definitions in order to identify diffractive events. 
Although these effects are found to be significant, it is shown that once they are 
properly taken into account, a very precise agreement between diffractive DIS and 
diffractive dijet photoproduction emerges without any significant hint of hard fac- 
torization breaking. 
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Introduction 



Recent years have seen a renewed interest in diffractive physics and specif- 
ically, in the possibility of using hard diffractive processes in order to inves- 
tigate 'at parton level' diverse features of diffraction [p]]. 

The success of Regge theory in the phenomenological description of soft 
diffraction || has naturally encouraged most analyses to use not only the 
language but also many critical hypothesis of that approach but in the hard 
regime and including simultaneously elements of perturbative QCD|l], |3|, §]. 
In this way, issues like the 'parton content of the pomeron' or the 'factoriza- 
tion properties of the pomeron parton densities' have become the center of 
active experimental programs and intense theoretical debates 0. 

This concoction of perturbative QCD and Regge ideas has shown to be 
quite successful, from the phenomenological point of view, when applied in- 
dividually to an increasingly large set of diffractive deep inelastic scattering 
data || 0] , at least as a first approximation, and also to diffractive processes 
within hadron-hadron collisions ||. However, it seems to have serious prob- 
lems when it is meant as a processes independent description of both sets 
of data |J. This fact could indicate a breakdown either of universality in 
the ingredients coming form the Regge picture, due for instance to an unex- 



pected energy dependence in the pomeron probability density |TT| , or in the 
hard factorization property assumed in the QCD inspired component of the 
treatment @, 0, |. 

Alternatively to the above mentioned approach, many hard diffractive 
processes, including diffractive DIS, can be treated rigorously within pertur- 
bative QCD -and without caring about the validity of the hypothesis implicit 
in the Regge approach- using the framework of fracture functions [0. This 
framework was originally thought to deal with semi-inclusive processes within 
QCD in a fully consistent way fll3f , and more recently it has been shown to 
be particularly successful in the phenomenological description of diffractive 
DIS and leading baryon production, and also in the analysis of departures 
from the standard Regge hypothesis, already found in the available data [|TJ] . 

In the present paper, the framework of fracture functions is applied to 
diffractive photoproduction of dijets |TJ], This is meant not only as 
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a test of the general approach and as a complementary verification of the 
parameterizations extracted from DIS, but also as a test of the hard factor- 
ization property assumed for fracture functions in the case of hadron-hadron 



collisions. At variance with the case of diffractive DIS ||16|| , the factoriza- 
tion property of fracture functions, which is crucial in order to compute the 
sizeable contribution of resolved photons to the photoproduction processes 



14, [Lql, may be challenged by the occurrence of soft interactions before 



the hard scattering takes place in diffractive hadron-hadron collisions || [K]]. 

Performing this analysis, it has been found that the use of non-equivalent 
rapidity gap definitions in order to identify diffractive processes in the differ- 
ent experiences conspires strongly against the above mentioned tests. In the 
language of fracture functions, to adopt different gap definitions implies that 
only certain components of the whole fracture function are tested, suppress- 
ing or enhancing contributions coming from certain kinematical regions, and 
leading seemingly to a factorization breakdown. However, here it is shown 
that these effects can be properly taken into account yielding a very pre- 
cise agreement between diffractive DIS and diffractive dijet photoproduction, 
without any significant hint of hard factorization breaking. 

In the following section, the standard procedure used in the extraction of 
diffractive parton densities or, alternatively, fracture functions in diffractive 
DIS is reviewed giving emphasis to the issue of the use of different defi- 
nitions for the diffractive event. Then, we proceed to the computation of 
diffractive photoproduction of dijets analyzing the consequences of different 
conventions about rapidity gaps for this kind of process. Finally, results from 
the computation are compared with available data from the ZEUS and HI 



collaborations [O, ITH] and conclusions are presented. 



Parameterizations: 



It is customary to analyze diffractive DIS in terms of a so called 'diffractive 
structure function' F^ 3 \ defined through the triple-differential scattering 
cross section 
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in the usual kinematical variables 
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where q, P, and P' are the momenta of the virtual boson, of the incident 
proton, and of the final state proton, respectively, and y = Q 2 /(xs). Both 
the cross section and the structure function defined in this way, include an 
implicit integration over a given range of the variable t = (P — P') 2 , and in 
order to identify diffractive events, different selection criteria are commonly 
used. From a theoretical point of view, the most straightforward method is 
the positive identification of a final state proton carrying a substantial frac- 
tion of the incoming proton momentum. However, and mainly for technical 
reasons, the most widely used method retains events with a pseudo-rapidity 
gap (absence of hadronic activity) between the two hadronic systems in which 
the hadronic final state is conventionally divided ||. 

In the Regge inspired approach, the diffractive structure function is as- 
sumed to be given by the product of the probability f F / p (x ]P ,t) to find a 
pomeron in the incoming proton, which only depends on the variables x F and 
t, and a 'pomeron structure function' F^(/3,Q 2 ), given by parton densities 
which behave according to Altarelli-Parisi evolution equations and factorize 
as ordinary parton distributions M. Aside from the issue of the validity of 
these assumptions, it is clear from the preceding paragraph that non equiva- 
lent selection criteria for diffractive events can lead to the extraction of non 
equivalent pomeron probability densities and pomeron parton distributions, 
eventually non universal |jl8| . 

Alternatively, in reference it has been shown that the diffractive 
structure function defined as in Eq.(l) just represents the low x p limit of 
a more general semi-inclusive process, which in perturbative QCD can be 
rigorous and thoroughly described using fracture functions. 

In leading order, fracture functions account for target fragmentation com- 
ponents in the full semi-inclusive deep inelastic lepton-proton scattering cross 
section with an identified proton in the final state. These components are 
neglected in the most usual approach based only in ordinary proton parton 
distributions qf (x, Q 2 ) and proton fragmentation functions Df (z, Q 2 ), which 
accounts for current fragmentation processes. 
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Then, the full semi-inclusive cross section, 

d 3 a p/p Ana 2 




i-y + ir x (3) 



dx dQ 2 dz x Q 4 
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includes an additional contribution that can be written as a sum of hard 
scattering cross sections weighted by the fracture function densities corre- 
sponding to the different quark flavors i for protons fragmenting into protons 
x, z,Q 2 ), which can then be understood as the probabilities to find a 
parton of flavor % in an already fragmented proton [pL3| . 

In Eq.(3), the variable z is, as usual, the ratio between the energy of 
the final state proton and that of the proton beam in the center of mass 
of the virtual photon-proton system. For very forward protons z ~ 1 — x p 
and the current fragmentation components in the cross section fade away 
highlighting the relation between fracture functions and diffractive or leading 
proton structure functions. 

As for ordinary structure functions, although the scale dependence of 
fracture functions can be computed within perturbative QCD, they are es- 
sentially nonperturbative objects which can not be obtained from first prin- 
ciples in a perturbative analysis. However they can be extracted from global 
QCD analysis of experimental data. In reference |12] it has been shown that a 



very simple parameterization in the variables (3 and x p and at an initial scale 
Ql = 2.5 GeV 2 for singlet quarks and gluons in the corresponding fracture 
function, 

xM? /p ( Xjp ,P,Ql) = NiF^l-P)* x (4) 

{/3a£ + C LP {1 - /3)^[1 + a LP (l - x p )^]} 

reproduces both diffractive and leading proton data as complementary regimes 
of the same semi-inclusive process without need to rely in Regge model as- 
sumptions. There, the values for the parameters were obtained fitting simul- 
taneously diffractive and also leading proton positron-proton DIS data from 



the HI Collaboration || [19| and can be found in reference [fT2 |. The former 
data fix the low x p behavior of the fracture function while the latter do the 
same but for larger values of x p . 
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Regarding the diffractive selection criterion for DIS data, HI retains 
events having a rapidity gap that spans at least the region 3.3 < r\ < 7.5, 
assuming that this gap definition guarantees the dominance of single diffrac- 
tive (semi-inclusive) events like e + (k) + p(P) — > e + (k') + X + p(P'), where 
X is a generic hadronic final state, and p(P') the final state proton |J. This 
assumption in fact seems to be justified, at least for this kind of process and 
in the kinematical regime of the experiment, given the remarkable agreement 
between the outcome of the parameterization and data obtained by ZEUS 
[0] detecting the final state proton, as pointed in [ 12| . However, the assump- 
tion not necessarily holds for other gap definitions, processes, or kinematical 
regimes. In the following section we will show that diffractive dijet photo- 
production provides an illustrative example of this situation. 



Photoproduction 



Both the ZEUS and HI Collaborations have presented results of diffrac- 



tive dijet photoproduction experiments at HERAjlv], [14|, [Tq| . These pro- 
cesses, also depicted by e + (k) +p(P) — > e + {k') +X +p(P'), are characterized 
by the presence of at least two jets in the hadronic final state X that ac- 
companies the outgoing proton p(P') and positron e + (k'). The diffractive 
criterion is implemented imposing upper limits for the pseudo-rapidity rf^ x 
of the most forward particle belonging to the hadronic state X and with 
energy in excess of 400 MeV {v^ x < 3.2 in the case of HI, while ZEUS takes 

'/,';;;;', < i-s). 

The photon interaction in these processes can be either of direct or of re- 
solved nature, making convenient the introduction of the variable x 7 , which 
measures the fraction of the photon momentum that takes part in the hard 
interaction. The direct component of the diffractive dijet cross section can 
be straightforwardly computed replacing in the inclusive (non-diffractive) 
cross section, the usual proton parton distributions by the corresponding 
fracture densities xM l p , (x F , (3, /i 2 ) and allowing an integration over the per- 
tinent range in x p pU[ , 
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where f 7 / e (y) is the flux of photons from the positron, px the transverse 
momentum of the outgoing partons, /i the scale at which the strong coupling 
constant is evaluated, and a i+1 ^k+ P the parton-photon cross section with two 
partons in the final state. 

For the resolved component, the procedure would be analogous to the 
preceding one with the proviso that in this case no formal proof of hard 
factorization guarantees it, 



a 



resolved _ / dy f ^ j j ^ j d p xM ^ j rfp ^+g"^ . (6) 



Here . J~ j j ry sy j i^l 2 ) denotes parton distributions in the photon, <j i+ j^ k+p the 
parton-parton cross section with two partons in the final state. For simplicity 
we have dropped the arguments for the densities and also the sums over the 
partonic indexes i and j. 

In any case, before dealing with the issue of hard factorization and making 
comparisons between the parameterization and the available data sets, it is 
imperative to analyze the compatibility between the gap criteria implicit in 
them. 

Given that the profile of the hadron activity as a function of the pseudo- 
rapidity of the most forward hadron, and thus the consequences of a given 
gap convention, depends on various non perturbative hadronization effects, in 
what follows we simulate them using a variant of the Monte Carlo generator 



POMPYT 2.6 [21 1, modified in order to include parameterizations for fracture 
functions. 

The inclusion of fracture functions in POMPYT imply to abandon explic- 
itly the pomeron flux factorization hypothesis, modifying both the shape and 
scale dependence of the diffractive parton densities as a function of x p , as 
required by DIS data [Q. In what follows, all the simulations are computed 
using the best fit to diffractive DIS data, denoted as SET A in reference |L2] 



as the fracture function, and SaS 2M from reference |22 for the photon. 



In Figure (1), dijet photoproduction events generated in the kinematical 
regimes of HI (Fig. la) and ZEUS (Fig. lb) are shown as a function of 
Vmax- The thick vertical lines represent the upper cuts in rj^ x applied by 
each collaboration. 
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Figure 1: Dijet photoproduction events generated in the kinematical regimes 
of the a) HI and b) ZEUS experiments 



It is worthwhile noticing that even for the same kind of process, the 
different kinematical regimes covered by both experiments given for example 
by cuts in the transverse energy Et and pseudo-rapidity i] of the two most 
energetic jets, etc., yield rather different distributions in rft^x an d thus, are 
affected in a different way even for the same pseudo-rapidity gap definition. 
For ZEUS, approximately 25 % of the total number of events are concentrated 
at rftntx < 1-8, while for HI kinematics approximately 10 % of them survives 
the same constraint. 

In the computation, an additional multiplicative normalization factor has 
been applied to the fracture function in order to take into account the fact 
that the original parameterization was adjusted to data already filtered by 
a gap definition. This normalization factor can be defined in a first ap- 
proximation as the ratio between the total number of dijet photoproduction 
events and those that satisfy the gap definition implicit in diffractive DIS 
data, both quantities computed in the kinematical range where the simula- 
tion takes place. In the present case, this rate is found to be 1.14 for ZEUS 
and 1.61 for HI. 

Figures (2) and (3) show the comparison between dijet photoproduc- 
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Figure 2: Comparison between HI dijet photoproduction data and the outcome 
the fracture function parameterization. 



tion data [TJ[ and the outcome of the fracture function parameterization 
obtained from DIS. The cross sections are computed in picobarns and as 
functions of the hadronic level variables defined within each experiment. 

As it can be seen, the agreement is quite good for both sets of data 
although the effects of taking into account the different rapidity gap con- 
ventions are large. It is important to notice that neglecting the stringent 
restriction applied in the ZEUS measurement {rft^ x < 1.8), the computation 
would have led to cross sections up to five times larger. 
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Figure 3: Comparison between ZEUS dijet photoproduction data and the outcome 
the fracture function parameterization. 

The comparison with ZEUS data shows that the numerical results tend 
to underestimate the cross sections for large values of f3 obs . This fact could 
indicate the need of a further adjustment in the diffractive parameterization 
in that region, which in fact is allowed by the present level of uncertainties 
in DIS data, specially at lower scales. 

For HI photoproduction data, the gap definition is very similar to that 
of diffractive DIS data so the effect of neglecting at all such cuts is rela- 
tively small when analyzing event distributions in variables other than f]^ x . 
However, if one wants to reproduce the whole picture including the rf r 
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dependence, of course they can not be neglected and some care must be 
taken in order to not apply the cuts twice, as pointed by the large normal- 
ization factor. 

In the precedent discussion, specifically in the computation of the nor- 
malization factor, we have assumed that the rapidity gap restriction only 
reduces uniformly the number of events, i.e., without changing the distri- 
bution of events in other kinematical variables. This seems to be a good 
approximation for the Et and rj distributions, but it is not for x p . The 
Monte Carlo simulation shows a clear correlation between the cuts applied 
in r)^ x and the mean for the distributions in x p implying that the larger is 
the gap, the lower is the mean x p proved, as shown in figure (4). 
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Figure 4: Distribution of dijet events as a function x F for different rapidity 
gap definitions (a) < 1.0, b) rf^ x < 2.0, and c) < 3.0) 



The direct consequence of the above mentioned effect is that different 
gap definitions weight differently contributions coming from a given range of 
x p , thus adding to the discrepancy between diffractive structure functions 
extracted from different experiments. In our approach, this inconvenience 
has been straightforwardly overcame computing the normalization factor for 
various bins in x p . The main effect of this procedure in the comparison with 
HI data is to enhance the cross section at low z p (large x p ) in figure (3d) 
and reduce it slightly at large z p . For ZEUS, the normalization effect is small 
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so the correction is almost negligible. 



This considerations about rapidity gap definitions apply also to the com- 
parison between diffractive DIS and diffractive dijet production within proton- 
antiproton collisions at Tevatron, where large discrepancies have been ob- 



served H, [1T[] . In the comparison of DIS data and these processes, the effects 
associated with the gap definition can even be enhanced not only because of 
the differences in the kinematical ranges covered but also in the extrapolation 
of the diffractive DIS parameterizations to rather large values x p , dominant 
in proton-antiproton collisions. 



Hard Factorization 



Regarding the issue of hard factorization, it is worth noticing that in 
both ZEUS and HI photoproduction experiments, and even in the regions 
dominated by the contributions coming from resolved photons (x 7 < 0.7), 
no hints of significant factorization breaking effects are found, confirming 
our initial assumption about hard factorization in these kind of processes. 
This is particularly apparent for HI data which contains a much larger frac- 
tion of events initiated by resolved photons and would indicate that hard 
factorization, if not exact, is a very good approximation. 

It is also interesting to notice that at variance with the analysis of HI 
data of reference |15||, in the present treatment no substantial overestimate 



of the cross section as a function is found. In reference [151 this excess 



was interpreted of a reduced gap survival probability |23|| , and an ad hoc 



correction factor was fixed with the data. In the present approach these 
discrepancies can be traced back to the non equivalence of rapidity gap defi- 
nitions, which encompasses effects related to spectator interactions, but in a 
more predictable way. 



Conclusions: 



Dijet diffractive photoproduction has been analyzed in the framework 
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of fracture functions finding that the main effects that conspire against a 
unified picture for diffractive DIS and diffractive photoproduction in terms 
of diffractive parton distributions are related to non-equivalent criteria used 
to select diffractive events. 

The connection between different rapidity gap definitions is model-dependent 
in the sense that it includes non-perturbative hadronization details. How- 
ever, using a slight modification of the POMPYT Monte Carlo generator, 
devised in order to include in it fracture functions, a consistent and accurate 
picture can be drawn. 

In this way, the analysis suggests that the supposed discrepancies between 
the diffractive parameterizations coming from these experiments are rather 
due to the non-universal character of rapidity gaps, than to a hard factoriza- 
tion breakdown, and that the potential factorization breaking mechanisms 
are negligible, well beyond the accuracy of the present data. 
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